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■ Abstract. We consider a nonlinear model that is a combination of the anisotropic two- 

dimensional classical Heisenberg and Toda-like lattices. In the framework of the Hirota direct 
' approach, we present the bilinear equations describing this model and derive its N-soliton 

rS^ , solutions. 
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1 Introduction. 

In this paper we consider the two-dimensional lattice that can be viewed as a generalization of 
the anisotropic two-dimensional classical Heisenberg model [U [2] , 

Cn ; ^Heis = X] X] i^r' ^r+Si) ■ (1-1) 

^ ! reA«=l,2 

^SJ ' where (p^, is a three-dimensional unit vector, 

o ■ 

^: (0.,</'.) = i (1.2) 

(with brackets standing for the standard scalar product), Ji^2 are the constants characterizing 
^ I the interaction between near-neighbour sites (exchange constants) and A is a two-dimensional 

H ' lattice formed by two vectors Si and 

A = {mi^i + rn2^2}^,,^,=o,±i,±2,... (1-3) 

The generalization that we are going to study consists in replacing the constants Ji with some 
functions of new variables. In more details, we associate with each site, in addition to the vector 
4>r^ a new variable Ur and modify the exchange constants as 

Ji ^ JieyiY>{ur - Ur+Si) ■ (1-4) 

An elementary example that leads to the above modification of the exchange interaction is 
to permit the spins to oscillate in the direction perpendicular to the plane and to state that 
the interaction coefficients J depend on the distance (in the three-dimensional space) between 
the spins: £ab = J {\Ra — Rh\) {4>ra^4'rt) where Ra = ra + Ur^v, with i/_L5i^2- In this case 
the interaction between the nearest neighbours depends on \6i\ and | I, that can be 

modeled by (jl.4p . Of course, the dependence given by (|1.4p is far from being realistic, however 
this toy model can give some insight into effects caused by such kind on nonlinearities, and 
especially into the possibility of appearing of specific structures like solitons that are discussed 
in this paper. 
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Figure 1. Direct and dual lattices, nearest-neighbour notation. 



To summarize, our model is described by the energy functional 

f = ^ ^ Ji exp {Ur - Ur+S,} {(f>r, (t>r+d,) ■ 



(1.5) 



rGAi=l,2 



It is easy to see that neglecting the (p^-part, or imposing the restrictions (0^, 4'r+6i) — fo'^ 
r and redefining the constants Jj , one arrives at the one of the Hirota's versions of the discrete 
2D Toda lattice [31 S]: 



rgA i=l,2 



(1.6) 



whose field equations 6£Toda./Sur = are known to be integrable (see, e.g., [5]). Thus, we call 
model (II. 5p . which is the subject of this paper, the two-dimensional Toda-Heisenberg lattice 
(2DTHL). 

To make the following formulae more readable we introduce the alternative notation: instead 
of the vector index we will use a letter one, 



(1.7) 



and denote the nearest neighbours of the point A as indicated in figured] (R, L, U and D stand 
for 'right', 'left', 'up', 'down'). The energy of the model can be rewritten as 



where 



Sa = '^Jb exp {eb {ua - Ub)} {4>a^ 4>b) 

B 

and the summation index runs over the nearest neighbours, 

E-= E - 



B=R,U,L,D 



The constants Jb are nothing but Ji^2, 
Jr = Jl = Ji, Ju = Jd = -h 



(1.8) 



(1.9) 



(1.10) 



(1.11) 
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Table 1. Nearest-neighbour notation and parameters. 



while Eb take into account the signs of Ug in the exponential functions, 

En = -El = Eu = -Sd = I (1-12) 

(see table [T]). The central object of the presented study are the Euler-Lagrange equations 
corresponding to (jl.Sp with (jl.9p . 



T-^ = ^ ^bJb exp{eB {ua - Ub)} {4>ai ^b) = 0) (1-13) 



Sua 
5£ 



54> 



^Jb exp{eB {ua - Ub)} [0^ X (t)g] = 0. (1.14) 



A 



Here, we have calculated the derivative with respect to 0^ bearing in mind the fact that |0^| = 1, 
which implies that the admissible variations should be orthogonal to 0^, which in its turn leads 
to the 'definition' 5(</)^,i/)) /6cf)^ = [(|)^,^J;]. The same result, (jl.l4p . can be reproduced by 
means of the Lagrange multipliers. 



2 Bilinearization of the 2DTHL. 

One of the important moments of the direct approach is to introduce an appropriate represen- 
tation of the unknown functions, 0^ and Ua in our case, that leads to the bilinearization of 
the equations we want to solve. To this end we first pass from the vector variables 0^ to the 
complex ones, Qa and r^, defined by 



(2.1) 




(2.2) 



(2.3) 

(the asterisk denotes the complex conjugation). In terms of Qa and the scalar and vector 
products are given by 

(0^,0 J = — ^(1 + Reg^g*), (2.4) 
\/PaPb 

^ lm{qA-qB)\ -i \ (2.5) 
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and one can present the left-hand sides of our main equations, (I1.13P and (jl.l4p . as 

oua \'PaPb 



and 

SSa _ 1 
54)^ 2i 

where 



(2.7) 



:= EJs'-^^^^^^^^i,.-,.), (2.8) 

7^. := ^Js'-^^^^^p^^irA-r.). (2.9) 
Thus, the field equations of our model become 

VA = QA = nA = 0. (2.10) 

Considering the functions n^, or Ur, we first choose the 'boundary conditions', 

Ur ^ {k,r) as |r| — oo (2.11) 

where is a constant vector determining the asymptotics of solutions, and present them as 

Ur=Ur + {k,r) (2.12) 

with bounded Ur- This change of variables modifies the exponential functions that appear in 
^T^ . ^B) and (1231) as 

exp {eb (ua - Ub)} = exp {e^ (n^ - Ub)} (2.13) 

with 

^fl = = 6 = exp {- (fc, (5i)} , 

(2.14) 

Cu = Cd = ^2 = exp {- (fc, 62)} 
(see table d]). 

Now we arrive at the key moment of bilinearization of our equations. It consists of introducing 
the dual lattice, as is shown in figure [H whose nodes closest to the point A are denoted by the 
letters A^, W, S and E (coming from 'north', 'west', 'south', and 'east') and the tau-functions 
T, a and p. 

The real functions r are defined by 

Pa = "-^ (2.15) 

n 



and 



exp(2'u^) = — (2.16) 

Tn 
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while the complex ones, a and p are defined by 

= ^, = ^ (2.17) 

Ta Ta 

and are related by 

tI = ctaPa+TnTs (2.18) 

which is the consequence of the identity Pa = I — Qa^a- 

By simple algebra one can obtain the following representation of the exponential functions 
from ([2T3]) : 

exp {eb {ua - Ub)} = ^/PaPb ^"^^^ (2.19) 

where B' and B" are the ends of the oriented edge of the dual lattice that crosses the edge {AB) 
of the direct one (see table [T|). 
Equation (|2.19p leads to 

Va = ^e.J.e. ""^""^^^'^^"" , (2.20) 

Qa = VJ.^.^^, (2.21) 

^ Tb'-Tb" 

TZa = ^JbCb^-^^ (2.22) 
g Tb'Tb" 

where we use the shorthands 

{x,y}AB = XaVb + XbVa, [x,y]AB = XaVb - XbVa- (2.23) 
The final step in bilinearization of the field equations (j2.10p is the following ansatz: 

[^^t]ab = fJ-B[cr,T]j^,g„ (2.24) 
[p^t]ab = fJ-B[p,T]g,gn (2.25) 

with constant /ig. To make this ansatz consistent with the lattice translations one has to impose 
restrictions on these constants, namely 

p, = Pn = PU (2.26) 

that follow from the shifts along the direct lattice, and 

PiP2 = -1 (2.27) 

that follows from the shifts along the dual one. So, all the constants Pb can be written in terms 
of one of them, p = /i^j, as is shown in tabled! 

Ansatz ([TM]) . (p;^ transforms the quantities dMU, (j2:22|) into 

Qa = ^ JbCbPb {Qb' - Qb") , (2-28) 

B 

TIa = 5^ JflCs^fl(^fl' -^-B") (2.29) 
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and it is easy to see that to solve equations Qa = TIa = one has just to put 

Ja^BfJ-B = I* = constant. (2.30) 
Indeed, this equation, which can be rewritten as 



+ J2C2 = 



(2.31) 



and which determines as a function of Ji^2 and ^1^2 (or, equivalently, as a function of Ji 2 and 
k) leads to 



Qa = hY^giqB' - qB") =0, 
T^A = hY^girB' -rg") =0 

where the last identities follow from the cyclic character of B' and B" (see table [1]) . 
To solve the last field equation, Va = 0, one can use 



1 - ^{q,r}AB =Pa + 5/^, 



and rewrite Va as 

Va=V^a^ +V'a^ 



TaTb 



{q,r}AB' - {qir}A 



where 



{q,r}AB' - {q,r}. 



Expressing Jb£,b from (j2.30p one can get 



= (T^ - 1) J— [^,TuTs + ^1~\LT^) 

TaTw 

where indicates the shift in the 'east' direction: 

Tb/a = /b, Te/w = /a, Te/l = fs, etc. 
(2) 

In a similar way, Va can be presented as 

V^^^/h = {TE-l){q,r}Aw 
which means that equation Va = can be replaced with 

— - — ({cr, p}aw + fJ'TuTs + fJ'^'^TLT,^) = A = constant. 
TaTw 

The simplest way to meet the condition of boundedness of Ur is to choose 
Thus, we have proved the following 



(2.32) 



(2.33) 

(2.34) 

(2.35) 
(2.36) 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



(2.41) 
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Proposition 2.1. A wide range of solutions for the field equations ^.13\) and ^1.14^ can be 

obtained from the bilinear system 

= tI- TnTs - (JaPa 
= [(T,t]al + ^J'W,r]sw 

(2.42) 

= [p,t]al + l^[p,r]sw 

= {a, p}aw + fJ'TuTs + P~^TlTm - (/U + /U~^) TaTw 

This system is the one we were looking for: a bilinear system providing solutions for the 
field equations (jl.l3p and (jl.l4p . However, these equations, except the first one, are four- and 
five-term ones, which makes their solution rather cumbersome. It is possible to prove directly 
that solutions presented below (see section U]) satisfy (j2.42p . However we take another way and 
split, in the next section, this system in a set of standard Hirota-like three-term equations. Of 
course, such splitting narrows the class of solutions. Nevertheless, this class is rather rich and 
contains the soliton solutions that we want to obtain. 



3 Reduction of (12.421) to Hirota-like equations. 



To explain the idea behind the splitting we are going to do, it seems reasonable to pass from 
the vertex notation to one based on the shifts/translations T, 

/a = /, fw = '^wf, fs=Tsf, etc. (3.1) 

In this notation, the bilinear system (j2.42p can be rewritten as 

' = T^-pa- (Tsr) (T^r) 

= r (T^a) - (J (T^r) + /X pgr) (T^-(j) - (Tsd) (Th^t)] 

= r(T^p)-p(T^r)+/i[(Tsr)(T^-p)-(Tsp)(TH,r)] 

= a (T^.-p) + p (Th.(j) - (/i + p-i) r (T^t) + p (T^r) (T^r) + p-^ (T^r) (T^r) 

Each point of both direct and dual lattices can be reached by a composition of two basic shits, 
say, Ts and T^i,-, 

T^ = T^\ Ti = TsTvy, Tu = T~^Tw, etc. (3.3) 

and hence all equations can be presented in terms of these two translations only. The trick that 
leads to the reduction to three-term equations is to use three shifts as a basic system. One of 
them is, say, Tg while two more shifts come from the splitting of T^, 

= T^Ty (3.4) 

or, alternatively, T^Tv^- = Tx^y- Application of this construction to our equations leads to 

Proposition 3.1. A wide range of solutions for the field equations il.l3\) and ( [i. j^[ ) can be 

obtained from the bilinear system of Hirota-like equations 

= - pa - (T^r) (T^r) (3.5) 

= TiTxT)-a{Txp)-iT,sT){T^xr) (3.6) 

= p{Tsa){Txp)-{TsT){TxT) + T{Tsxr) (3.7) 

= TiTyT)- p{Tya)-iTsT){T^yT) (3.8) 

= p{Tsp){Tya)-{TsT){TyT)+T{TsYT) (3.9) 

One can find a proof of this statement in appendix O Geometrically, this constructions can 
be interpreted as if we considered our plane as a part of three-dimensional lattice and presented 
our bilinear equations as a projection of more simple three-dimensional system. 
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4 Soliton solutions for (^^-(^M)- 

In this section we derive the sohton solutions for the bilinear system (|3.5p - (|3.9p . This will be 
done in two steps. First we solve it without taking into account the condition r = —q* and, 
secondly, find the restrictions that should be imposed on the parameters of the solutions to meet 
it. 

It should be noted that equations ()3.5p - (l3.9p are nothing but the Ablowitz-Ladik system 
[H [71 [8]. Thus, we use some of the results of papers [Gj [7|, namely the structure of the soliton 
solutions without developing the inverse scattering transform from scratch. 

The main blocks for constructing the A^-solitons are the matrices that satisfy the 'almost 
rank-one' conditions 

LA - AR = |l)(a|, RB - BL = |l)(b|. (4.1) 

Here L and R are constant diagonal matrices, |1) is constant A^-component column, (a| and (b| 
are A^-component rows depending on the coordinates of the problem. 

We look for solutions of the form, similar to the form of soliton solutions for the Ablowitz- 
Ladik model: 

r = det|l-AB| =det|l-BA| (4.2) 

and 

g= (a|FR~^|l), r = -(b|L"^G|l). (4.3) 

Here 1 is the N x N unit matrix, 

F = (1-BA)"\ G = (1-AB)-^ (4.4) 

and (a|M|l) is the standard row-matrix-column product: (a|M|&) = '^jf.CLjMjkbk- where aj, 
Mjk and bk are components of a row (a|, a matrix M and a column 

The second part of the 'solitonic ansatz' is that the action of the shifts Tg^x,Y,... can be 
implemented as the right multiplication by diagonal constant matrices. The structure of these 
matrices, that are rational functions of L and R, can be obtained, again, from [6l [7]. The 
resulting formulae can be written as follows: 



Proposition 4.1. Soliton solutions for 11 3. 5\) - [379]) are given by 

r = det|l-AB| =det|l-BA| (4.5) 

and 

a = r(a|FR~^|l), p = -T{b\l"^G\l) (4.6) 
with the following actions of the T -shifts: 

TgA = AR TgB = BL"^ 

T;,A = A{l+fiR) Tj,B = B(l+/iL)-i (4.7) 
TyA = A(l+;uR-i)"^ TyB = B(l+^L-i). 
Equations (14. 7p lead to similar ones for the rows (a| and (b|: 
Ts(a| = (a|R T,.(b| = (b|L-i 

T^-(a| = (a|(l+/iR) T^(b| = {b\ {1 + i^iL)-^ (4.8) 
Ty(a| = (a|(l+^R-i)~^ Ty(b| = (b| (l + ^uL^i) . 
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This proposition is proved in appendix IbI 

After having derived the 'general' soliton solutions for (I3.5p - ()3.9p . we have to ensure the 
'physical' involution r = —q*, or 

r* = r, p = -a* (4.9) 

bearing in mind that 

^* = e/i, e = ±1 (4.10) 

which follows from the fact that = —J2(,2/JiCi is real, but possibly negative, number (see 
(12311) ). 

Omitting rather straightforward calculations, we present here the following results: the rela- 
tionships between the matrices A and B, and hence between the rows (a| and (b|, that ensure 
TM are 



L* = eR-\ R* = eL"^ (4.11) 
together with 

A* = -RBL-^ B* = -LAR-^ (4.12) 

and 

(a*| = e(b|L-2. (4.13) 
Restrictions (|4.12p - (j4.13p can be resolved by introducing, instead of A and B, one matrix C, 

A = CR, B = -eR-iC*, (4.14) 
and one row (c| instead of (a| and (b|, 

(a| = (c|R, (b| = (c*|L. (4.15) 
Definitions (j4.14p and (|4.15p enable to present F and G as 

F = R-iHR, G = H* (4.16) 

where 

H = (l + eC*C)-^ (4.17) 
and to rewrite (j4.2p and (j4.3p as 

r = det|l + eC*C| (4.18) 

and 

g=(c|H|l), r = -(c*|H*|l) (4.19) 
which clearly demonstrates the fulfilment of ()4.9p . 
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5 Solitons of the 2DTHL. 

Now we have all necessary to present the A^-soliton solutions for equations (|1.13|) . (jl.l4p . i.e. to 
write down the expressions describing the solitons of the 2DTHL. 

First we have to return from the next-neighbour notation to the 'absolute' one noting that 
for any lattice vector, 

r = midi + 171262, (5.1) 
and any function fr 

= T-iT™2/ (5.2) 
where / is the value of the function at some fixed point. Recalling that = T^Ty and that 

Tu = T~\ Tu = T^Tf (5.3) 
one can introduce the matrices C,. and the rows (c,.| by 

C^ = CM^lM^^ (c^l = (c| Mf^M^^ (5.4) 

where, as follows from (14. 7p . 

Ml = (l + /xR-i)(l+/^R)-\ (5.5) 
M2 = (R~i + ^a) (R + ^1)^^ (5.6) 

with constant C and (c| related by 

LC-CR = |l)(c| (5.7) 

or 

(C). _ Cfc _ eR*Ck ^^^^^ 
^ Lj — Rk 1 — eRjRk 

The parameter fi depends on vector k that defines the asymptotics of solutions (see (|2.3ip and 



A. = Mfc) = A/-^exp[i(fc,(5i-52)] (5.9) 



while the constant e is given by 

e = — sign J1J2. (5.10) 
In this notation, the main formulae of the previous section become 

g^ = (c^|H^|l) (5.11) 

where 

= (1 +eC;Cr)"^ (5.12) 

and 

= det|l + eC;C^| . (5.13) 
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Finally, one can present the soliton solutions of the model considered in this paper as follows: 
the vectors (pr are given by 

sin Or cos ifr 

Sm.OrSViHpr I (5-14) 



with 



cos Or 



6*,. = arcsin — ^, Lpr = B^igqr (5.15) 

1 + l^rl^ 

where qr are defined by (jS.lip while 

1 r+ 

= (fe,r) + -In^ (5.16) 

2 Tr 

with 

r± = T^Vr = det |l + eC;L^^C^R^M . (5.17) 



6 Conclusion. 



In this paper, we have presented a nonlinear 2D lattice and have obtained its A^-soliton solutions. 
Existence of such solutions means that the 2DTHL passes the 'three-soliton' test |10l [TH 1121 
[T3t [T5] (see also [161 [T7] for its application to discrete systems) which is a strong evidence 
of the integrability of the problem. However, the situations with the 2DTHL can be more 
complicated: there is a possibility that equations from Propositions 12.11 and 13.11 belong to the 
family of the so-called 'conditionally integrable' systems introduced by Dorizzi et al |18j . In any 
case, the questions related to the integrability of the 2DTHL, its zero-curvature representation 
and conserved quantities surely deserve further studies. Another possible continuation of this 
work is related to prospective applications to the theory of magnetic system. Traditionally, the 
soliton solutions are most representative ones for any integrable system. However, from the 
viewpoint of the theory of magnetism, it would be interesting to find other families of solutions, 
except the non-topological solitons presented in this paper, that correspond to more typical 
magnetic structures as, e.g., domain walls and finite magnetic domains. 
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A Proof of Proposition 13.1 



The proof of Proposition 13.11 is straightforward and consists in presenting the right-hand sides 
of the equations from Proposition 12.11 as combinations of the right-hand sides of the equations 
from Proposition 13. li To this end, consider the following definitions: 



a 

(,(1) 

r(2) 



T^-pa- (Tst) (T^r) , 

r {Tj,t) - a {Tj,p) - (Tgr) {T^j,t) , 

p (Tsa) (T^p) - (Tst) (T^r) + r (Tsxt) , 

T (TyT) - p {Tya) - (Tst) (T^yr) , 

p (Tsp) (Tya) - (Tst) (TyT) + r (Ts^r) 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 
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(which are shortcuts for the right-hand sides of (|3.5p - (j3.9p ). Our aim is to prove that vanishing 
of 0, b^^'^-* and c^^'^^ implies vanishing of 

n = T(Tr.a)-a{Tr.r)+fiiTsT)(T^a)-fi{Tsa){T^T), (A.6) 

^R = TiTr.p)- p{Tr.T) + fiiTsT){T^p)- fi{Tsp){T^-T), (A.7) 

= a{Twp) + p{Twa)- {p + p-^)t{Twt) (A.8) 
+p (Tst) (Tut) + /i-i (T^r) (T^r) 

(the right-hand sides of equations (j3.2p ). The first part of the proof is simple: it can be shown 
that and $H are linear combinations of b^^''^^ and c^^'^^. Indeed, one can check that 



{T^p)Q = (T,r)b«-(T^-r)b(2)-r(T^c«j+(T,r)(T^;,c(2)j^ 

{Tya)^R = -T(T,-bW)+(TsT)(T^,-b(2)) + (T,T)c«-(T^r)c(2). (A.IO) 



Thus 

^(1,2) ^ ^(1,2) ^ Q ^ £j^3^^0. (A.ll) 

As to *p, the calculations are slightly more complicated and can be performed in two steps. 
First, writing the system composed of (|A.2p and shifted (jA.Sp . 



b(i) = T{T^T)-a{T^p)-{TsT){T^^T) 

T^xc(2) = p (T^a) (Tx-p) - (Tw-r) {Tj,t) + (T^r) (T^^r) , 



(A.12) 



and eliminating (T^p) one arrives at 

CT (t^xc(2)) + p (T^-a) b(^) = A (T^r) - B (T^^t) (A.13) 

where 

A = /iT (T,^(j) - a (T„,t) , (A.14) 

B = // (Tst) (T^a) - a (T^r) . (A.15) 

In a similar way, elimination of Tj^xP from 

T^b(2) = paiTr,xP)-riT^xr) + {Tr,T){TxT) 
T^^c(i) = (T,,.t) (T^xt) - (T^a) (T^^p) - (T^r) (T^r) 



(A.16) 



leads to 

(T^a) (T^b(2)) + (t^;,c(i)) = -C (Txt) + (T^^^r) (A.17) 

with 

C = pa {Tut) - (Tjvt) (T^a) , (A.18) 

D = pa {T^-t) - T (T^-a) (A. 19) 



or, in the matrix form. 



A -B\ fTxT\_ fl,W 
-C D j It^^t j - U(2) 



(A.20) 
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where t)(^'^) are linear combinations of shifted b*^^'^^ and c^^'^\ On the other hand, one can 
straightforwardly verify the identity 

a {T^a) ^ + ^x-^{AD - BC) + {T^af a + (T,^o) = 0. (A.21) 

Equations <\2>.2\ imply vanishing of the last two terms of the left-hand side of the last equation 
and of the vector (t)^^\ t)^^^)"^, which means that the determinant of the matrix that appears in 
the left-hand side of (|A.20j) is zero, AD — BC = 0, since we assume r 7^ (and, hence, T^t 7^ 
and TjvxT" / 0). This, together with (lA.2ip . leads to 

a = b(^'2) = c(^'2) =0 ^ = (A.22) 



which completes the proof of the fact that equations (13.5p - p.9p imply (|3.2p which, in its turn, 
completes the proof of Proposition 13.11 

B Verification of the solitonic ansatz. 

To make the following formulae more readable we use here the 'node' notation not only for the 
nodes of lattices (direct and dual) but for the results of the auxiliary shifts Tx,y as well: 

fx=Txf, fy = Tyf, fsx = Txfs = TxTsf, etc. (B.l) 



First let us prove the fact that equations (I4.3p - (j4.8p imply (|3.6p . Applying (|4.ip and the 
S'-part of ()4.7p to the product BgAg one can get 

BA- BsAs = Bs|l)(a| (B.2) 

which leads to 

FF-i = 1 + FB5|l)(a| (B.3) 

and 

FsF-^ = l-FsBs|l)(a|. (B.4) 
Shifting this equation in the TjvTx-direction one can obtain 

KF^F-iK"i = l-KF^B;,|l)(a| (B.5) 

where 

K = R-i + ^ (B.6) 
(we have used the fact that (ajvx|K^^ = (a|). Multiplying ()B.3P and (jB.SP one arrives at 

FF-i KF,^F-i,K-i = l + FU|l)(a| (B.7) 

where 

U = Bs-F-^KF^B;, (B.8) 

= B<j(l+AsKF;,Bx)-KFxBx. (B.9) 

Noting that A5K = Ax and applying the identities 

FB = BG, 1+AFB = G (B.IO) 
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one can continue the calculations as follows: 

U = B,,(1 + A.^F^B^)-KF;,B^ (B.U) 

= (Bs-KB,OGx (B.12) 

= R-i|l)(b,.^|G;,. (B.13) 

After substitution of U in (|B.7p and calculating the determinants of the both sides one arrives 
at 

= l-qrx (B.14) 

(here, the identities r = 1/ det F and det(l + |u)(v|) = 1 + (v|u) have been used) which at the 
level of the tau-functions is 

= TTx - apx - TsTnx- (B.15) 



Thus we have proved that our solutions satisfy equations (|3.6p . 

Noting that the implementation of the shift Tx depends analytically on n and Tx|^=o 

the unit operator, one obtains from the above calculations, by sending ^ — t- 0, that our ansatz 
ensures (|3.5p as well. 

In a similar way, applying (14. ip to B^Ay one can obtain 

Fx^-^ = l-fifxBx\l){B\ (B.16) 

FF"^ = 1 + fiFBx\l){a\ (B.17) 

and then 

Fsfsx R~'F;,F-iR = 1 + /iF,.V|l)(a5| (B.18) 



where 



V = B^;, - F-^R-^F^B.y (B.19) 

= Bsx{l+^xfxBx)-R-^fxBx (B.20) 

= (Bsx-R-'B^)Gx (B.21) 

= R-i|l)(bs^|Gx (B.22) 



which leads to 

TTsx 
TsTx 

and 



1 - iiqsrx (B.23) 



= ficTsPx - rsTx + TTsx- (B.24) 

This proves the fact that our solutions satisfy (13. 7p . 

Considering the Ty-equations, we do not present here the calculations similar to ones dis- 
cussed above, leaving the verification of (13. Sp and (13. 9p to the reader. 



